Imaging electronic quantum motion with light 
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Imaging the quantum motion of electrons not only in real-time, but 
also in real-space is essential to understand for example bond 
breaking and formation in molecules, and charge migration in 
peptides and biological systems. Time-resolved imaging interro- 
gates the unfolding electronic motion in such systems. We find 
that scattering patterns, obtained by X-ray time-resolved imag- 
ing from an electronic wavepacket, encode spatial and temporal 
correlations that deviate substantially from the common notion 
of the instantaneous electronic density as the key quantity be- 
ing probed. Surprisingly, the patterns provide an unusually visual 
manifestation of the quantum nature of light. This quantum nature 
becomes central only for non-stationary electronic states and has 
profound consequences for time-resolved imaging. 

X-ray imaging | quantum electrodynamics | electronic wavepacket 

Abbreviations: TRI, time-resolved imaging; QED, quantum electrodynamics; 
DSP, differential scattering probability 

The scattering of light from matter is a fundamental phenomenon 
that is widely applied to gain insight about the structure of materials, 
biomolecules and nanostructures. The wavelength of X-rays is of the 
order of atomic distances in liquids and solids, which makes X-rays 
a very convenient probe for obtaining real-space, atomic-scale struc- 
tural information of complex materials, ranging from molecules Q] 
to proteins (2) and viruses (3). The power of X-ray scattering relies 
as well on the fact that the X-rays interact very weakly with the elec- 
trons in matter. In a given macroscopic sample, generally no more 
than one scattering event per X-ray photon takes place, the probabil- 
ity for multiple scattering being extremely small. The key quantity 
in X-ray scattering is the differential scattering probability (DSP), 
which is related to the Fourier transform of the electronic density 
p(x) as follows 0|5) 
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where is the differential scattering probability from a free elec- 
tron and Q is proportional to the momentum transfer of the scattered 
light. Procedures exist, for both crystalline (6) and non-crystalline 
samples Q, to reconstruct p(x) from the scattering pattern. Equa- 
tion <[X]) can be obtained from a purely classical description of elec- 
tromagnetic radiation scattered by a stationary electron density 0, 
yielding a result identical to that obtained from a quantum electrody- 
namics (QED) description of light. 

Equation Q} gives us access to a static view of the electronic den- 
sity. On the other hand, much progress has been made in recent years 
towards understanding electronic motion with time-domain table-top 
experiments lf8l[9] |10lllllll2|[T3l . owing to the availability of laser 
pulses on the sub-fs time scale 1 141 1 151 - An ultimate goal of imag- 
ing applications encompasses unraveling the motion of electrons and 
atoms with spatial and temporal resolutions of order 1 A and 1 fs, 
respectively |16|. Recent breakthroughs make it possible to generate 
hard X-ray pulses of a few fs |171|18]|, and pulses of length 100 as can 
in principle be realized fl9j [201 . Hence, a fundamental question that 
needs to be addressed is: How does an ultrashort light pulse interact 
and scatter from a non-stationary quantum system? 

Figure 1 shows schematically how the real-time, real-space dy- 
namics of an electronic wave-packet can be probed by scattering a 



short X-ray pulse from it. In the example selected here, the electronic 
wavepacket is a coherent superposition of eigenstates of atomic hy- 
drogen, with a wavepacket oscillation period of T = 6.25 fs. The 
electronic dynamics of the wavepacket is probed as a function of 
pump-probe delay time t. Previously, a similar example was used 
to illustrate time resolved imaging of an electronic wavepacket using 
ultrafast electron diffraction [211 . There, the instantaneous electron 
density was assumed to be the main quantity being imaged. 

In order to describe TRI, it is tempting to apply a similar treat- 
ment as used to obtain Eq. (TJ, i.e., a classical description of the 
X-ray pulse being scattered from a non-stationary electronic system 
described by the wavepacket. Therefore, under the assumption that 
the X-ray pulse duration is shorter than the time scale on which the 
motion of the electronic wavepacket unfolds (see supplementary in- 
formation), the expression for the DSP becomes 
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According to Eq. l[2}, TRI would be expected to provide access 
to the instantaneous electron density, p(x,t), as a function of the 
delay-time t. The idea to assign an additional degree of freedom 
to the electronic density in Eq. Q}, by replacing p(x) with p(x, t), 
is not new, and it has already been proposed and used in the recent 
past IU6H22ll23ll24l . However, a more careful consideration of the 
physics behind a time-resolved scattering process, quickly points to 
pitfalls in the simple logic underlying Eq. l(2j. Equation l(2j im- 
plies that the scattering of light does not affect the properties of the 
wavepacket. But in order to image quantum motion on an ultrafast 
time scale, one needs an ultrashort pulse, which has an unavoidable 
spectral bandwidth. Due to the finite bandwidth of an ultrashort 
pulse, it is energetically impossible to detect whether or not in the 
scattering process a transition between the eigenstates spanning the 
wavepacket, or transitions to other states closer in energy than the 
pulse bandwidth, has taken place. As a result, a physically correct 
treatment of ultrafast imaging from the electronic wavepacket will 
necessarily allow for transitions among electronic states caused by 
the scattering process (Compton type process). In other words, it 
must be expected that light scattering changes the wavepacket being 
imaged, and that this effect will be reflected in the image. In light of 
this, it is questionable whether Eq. ^ is at all justified. 

In the following, we apply a consistent QED description of light 
and a quantum mechanical description of matter. Only in this way 
can electronic transitions during the scattering event be properly ac- 
counted for 12511261 . We find that the photon concept, i.e. a quantum 
of light scattering from the electronic system, is crucial in the de- 
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scription of TRI. For a sufficiently short, coherent pulse, the resulting 
expression for the DSP is (see supplementary information) 

B-S/ A / W (»M°* , W> 1 ^. 
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Here, (••■}« is the expectation value with respect to the elec- 
tronic wavepacket at delay-time t. The operator C(H) = 
h J^f ln2 ex P(-aT^(H - (H)tf) is a function of the electronic 
Hamiltonian H and of the pulse duration r. AE refers to the pho- 
ton energy range accepted by the detector (see Fig. 2), and (H) t 
is the mean energy of the electronic wavepacket. ft(x) is the elec- 
tron density operator. The expectation value of n(x) with respect 
to the wavepacket at t gives p(x, t) i.e. (n(x)) t = p(x,t). As 
one can see from Eq. ((3), the quantum description of light pro- 
vides a very different expression for the TRI of the wavepacket, in 
comparison to that obtained from a classical description of light, 
Eq. l(2j- The crucial difference between both expressions is that 
Eq. (0 depends on (h(x.'))t(n(x)}t, whereas Eq. (0 depends on 
{n(x')C(//)n(x)) t . The fingerprint of electronic transitions within 
the finite spectral bandwidth of the X-ray pulse is encoded in C(H), 
which correlates the scattering of the photon at different times during 
the pulse duration. Similar interference effects caused by interaction 
of a wavepacket with a photon at different times during the pulse du- 
ration are known, for example, from fluorescence spectroscopy with 
optical lasers 1271 . Equations (2) and (3) rely on the assumption that 
the transverse and longitudinal coherence lengths of the incident X- 
ray radiation are larger than the size of the object. 

At this point it is interesting to note that many physical phenom- 
ena involving the interaction of light and matter can be explained 
from a classical description of light. Even the photo-electric ef- 
fect can be explained by assuming classical light 1281 , and very 
sophisticated experimental scenarios such as in photon antibunch- 
ing/anticorrelation experiments |29| 1301 or in Lamb shift measure- 
ments 1311 are necessary to prove the quantum nature of light 1321 . 
Surprisingly, a direct visualization of the quantum nature of light is 
obtained by looking at the time-resolved scattering pattern of a non- 
stationary quantum system. We illustrate the dramatic difference be- 
tween Eqs. (2) and (3) by presenting benchmark results for TRI of 
the electronic wavepacket introduced in Fig. 1. Due to the inher- 
ent finite bandwidth of an ultrashort pulse, it is necessary to compute 
the transition amplitudes from the electronic states involved in the 
wavepacket to all electronic states within the bandwidth. Depend- 
ing on the bandwidth, this can cover transitions to all bound states, 
including transitions to all Rydberg states, with any angular momen- 
tum. Therefore, the evaluation of Eq. <[3j even for hydrogen is nu- 
merically challenging. In contrast, the evaluation of Eq. I0 is easy 
as it involves only states within the wavepacket. 

Figure 3A shows scattering patterns, calculated with Eq. (3), in 
the Q x - Q z plane (Q y = 0) as a function of the delay time at times 
0, T/4, T/2, 3T/4, and T. Figure 3B shows the corresponding charge 
distribution of the wavepacket, which undergoes spatial oscillations 
along the z-axis. The probe-pulse duration in Fig. 3A is 1 fs. A 
photon energy detection width, AE, corresponding to 0.5 eV is used 
in the present calculation (see Fig. SI). Therefore, all transitions 
induced by the scattering process and resulting in a scattered photon 
energy within an energy range of 0.5 eV of the mean incoming photon 
energy were considered in the calculation of the scattering pattern. In 
order to ensure convergence, we considered transitions to states up to 
principal quantum number n — 40. Including all types of multipole 
transitions allowed by the conservation of angular momentum, this 
amounts to about 22000 states. 

As the electronic charge distribution oscillates towards the neg- 
ative z-axis, the pattern in Fig. 3A oscillates towards positive Q z 
values. At delay times T/4 and 3T/4, the electronic charge distribu- 
tions are identical, whereas the wavepacket carries a different phase, 



resulting in different patterns. It is evident from Fig. 3A that the scat- 
tering pattern is not at all times centrosymmetric i.e. it is not equal 
for Q and Q. This is in contrast with the centrosymmetric pattern 
expected from Eq. as a consequence of Friedel's law 1331 , and 
shown in Fig. 3C. The patterns shown in Fig. 3C were calculated 
using Eq. with a pulse of duration 1 fs. These patterns are more 
intense than the patterns shown in Fig. 3A, in which the intensity 
depends on the photon energy detection width. It is interesting to 
note that the pixel with maximum intensity varies as a function of 
the delay-time in Fig. 3A, whereas in Fig. 3C it remains fixed at 
Q = 0. This counterintuitive nature of the QED scattering pattern 
arises from the fact that it is not simply related to the Fourier trans- 
form of the instantaneous electronic density, but it is related to the 
electronic wavepacket through Eq. J3J. Clearly, the pattern in Fig. 3C 
is less rich than the correct pattern in Fig. 3A and it misses important 
phase information. Scattering patterns from electronic wavepackets 
corresponding to coherent superpositions of different sets of eigen- 
states (other than those used in Fig. 1) show similarly dramatic differ- 
ences between Eqs. ((2} and 0. Therefore, Fig. 3A is an unusually 
dramatic reflection of the quantum nature of light. 

In the illustrative example investigated here, an electronic 
wavepacket corresponding to the superposition of two hydrogenic 
eigenstates is considered. There have been several studies on the 
scattering of quasi-resonant light from strongly driven two-level sys- 
tems [34, 35 1, particularly under conditions when the system is pre- 
pared initially in a superposition of states 1361 . In those cases the 
transitions between electronic states are of dipolar nature. The TRI 
regime is, however, in the weak-field, one-photon limit and the pho- 
ton energy is highly detuned from the energy differences between 
states of the wavepacket. In TRI, scattering events are Compton-like, 
i.e., induced by the A 2 operator (see supplementary information), 
and can couple electronic states over a wide range of angular mo- 
menta. In order to illustrate the difference between these scenarios, 
we compute scattering patterns based on Eq. ([3} but restrict the calcu- 
lation to the two electronic states spanning the prepared wavepacket. 
The patterns in Fig. 4 present an intensity similar to that in Fig. 
3C because the same states are involved as in the calculation with 
Eq. 0. However, a careful look reveals the same periodicity and 
non-centrosymmetric features as seen in Fig. 3A. 

Coherent diffractive imaging (CDI) (37][38j, which provides ac- 
cess to structural information of periodic 1391 and non-periodic 1401 
samples, relies on the Fourier relationship between the scattering pat- 
tern and the electron density. CDI assumes Eq. for a station- 
ary electronic system, implying centrosymmetric patterns as a conse- 
quence of Friedel's law. Consequently, the natural extension of CDI 
to the time domain assumes that Eq. holds. On the other hand, 
we see that patterns obtained from a non-stationary quantum system 
are not centrosymmetric and cannot be related via Eq. to a real 
function p(x,f). Therefore, our results demonstrate that the exten- 
sion of CDI into the time-domain requires a careful analysis of the 
underlying processes. 

The scattering patterns calculated from Eqs. and are not 
only unrelated, but the effect of the pulse duration in the patterns 
is also totally different. The respective DSPs for a particular pixel 
in Q-space for different pulse durations are shown in Fig. 5. The 
curve calculated with the correct description, Eq. 0, has the same 
periodicity as the corresponding wavepacket dynamics (Fig. 5A). In 
contrast, the curve calculated with Eq. has the wrong periodicity 
(Fig. 5B). In case of pulses shorter than the characteristic time scale 
of the wavepacket, there is no optimal pulse duration for the signal 
computed with Eq. 0. The contrast as a function of time increases 
monotonically and for short enough pulses becomes almost constant. 
Conversely, the pattern computed with Eq. has an optimum con- 
trast as a function of time for a pulse length close to 1 fs. For shorter 
pulses, the pattern loses contrast again. The reason behind such an 
unexpected behavior is that in the case of hydrogen, the expectation 
value in Eq. ([3} becomes independent of t for r — > 0, intuitively be- 
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cause the electron has no time to change its position during r. Conse- 
quently, the pattern becomes independent of the electronic dynamics. 

We have shown that TRI from a non-stationary quantum system 
provides visual evidence of the quantum nature of light in a very 
simple way. Moreover, information on the quantum motion of an 
electronic wavepacket beyond the instantaneous electronic density is 
accessed by X-ray TRI with atomic-scale spatial and temporal reso- 
lution. Our result is counterintuitive as seen from the perspective of 
X-ray scattering from a stationary electronic system, as one would 
expect to access the instantaneous electronic density for sufficiently 
short pulses. Not only is this not the case, but the underlying physics 
fundamentally differs from such an interpretation. The illustrative 



example used here as a proof of principle lies in the time and energy 
range of interest corresponding to the dynamics of valence electrons 
in more complex molecular and biological systems 14 II 1421 |43l 1441 . 
We believe that our present findings will help to solidify the concep- 
tual foundations of the emerging field of TRI, where understanding 
the motion of electrons is key to understanding the functioning of 
complex molecular and biological systems. The advent of novel light 
sources will certainly bring us closer to this goal. 
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Fig. 1. Schematic of the time resolved X-ray imaging scenario used as an example throughout this work. An electronic wavepacket is prepared by a laser pump pulse 
(indicated in green) as a coherent superposition with equal population of the 3d and 4/ eigenstates of atomic hydrogen with projection of orbital angular momentum 
equal to zero. The polarization direction of the generated wavepacket is aligned with the laboratory z-axis. The electronic dynamics of the wavepacket is probed by 
an ultrashort X-ray pulse (indicated in blue) propagating along the y direction. A series of scattering patterns obtained by varying the pump-probe time-delay serve to 
image the electronic motion with high spatial and temporal resolution. 



Fig. 2. Schematic of the spectral energy distribution of Gaussian form of an incident ultrashort X-ray pulse with spectral bandwidth f u centered at the incident photon 
energy fUJk in ( see supplementary information). The window function WAE{h*>k s ) is a function of the scattered photon energy fiuJk s with a width AE (indicated 
by dashed blue lines). AE models the range of photon energies accepted by the detector. Therefore, elastic or inelastic transitions taking place during the scattering 
process that fall inside AE are considered in the evaluation of the scattering pattern (indicated by a green arrow). Transitions falling outside the range of AE are not 
considered in the evaluation of the scattering pattern (indicated by a red arrow). 
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Fig. 3. Scattering patterns in the Q x - Qz plane (Q y = 0) and electronic charge distributions of the wavepacket corresponding to the coherent superposition of 3d 
and 4/ eigenstates of atomic hydrogen (cf. Fig. 1). (A) Scattering patterns obtained with Eq. (5J, (B) electronic charge distributions and (C) scattering patterns 
obtained with Eq. [2], at pump-probe delay times 0, T/4, T/2, 3T/4, and T, where the oscillation period of the electronic wavepacket is T = 6.25 fs. The isosurface in (B) 
encloses ~ 26% of the total probability and has length 14-17 A along the z-axis and 7.5-9 A along the x and y-axes. The wavepacket is exposed to a 1 fs X-ray 
pulse with 4 keV photons. The patterns are calculated for Q max = 2 A -1 corresponding to 3.14 A spatial resolution. For the pulse used, a real-space pixel size of 
3.14 A requires the detection of photons scattered up to 60° with respect to the X-ray propagation direction, which is the j/-axis in the present case. The grid spacing 
in Q-space is 0.19 A -1 . The intensities of the patterns are shown in units of in both cases. 
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Fig. 4. Scattering patterns in the Q x - Qz plane (Q y = 0) obtained with Eq. [3] at pump-probe delay times 0, T/4, T/2, 3T/4, and T where the oscillation period of 
the electronic wavepacket is T = 6.25 fs. All the parameters used to obtain the patterns are the same as used in Fig. 3A. However, only the eigenstates belonging 
to the electronic wavepacket (3d and 4/ eigenstates of atomic hydrogen with angular momentum projection quantum number equal to zero) were considered. The 

d P 

intensities of the patterns are shown in units of -jrf. 
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Fig. 5. Effect of duration (full width at half maximum) of the X-ray pulse on the scattering patterns as a function of the pump-probe delay time. The electronic 
wavepacket oscillation period is T = 6.25 fs. The time evolution of the scattering pattern at an individual pixel as a function of the delay time is calculated using (A) 
Eq. (5J, and (B) Eq. (2), for different pulse durations, 100 as (black dots), 500 as (red dash-dots), and 1 fs (blue line). In the case of the pattern calculated using 
Eq. (3), there is an optimal pulse duration of about 1 fs, which provides a maximum contrast as a function of time. There is no such optimum for the pattern computed 
with the simple expression in Eq. |5J. 
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